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Abstract 

The Local Hamiltonian problem (finding the ground state energy of a quantum 
system) is known to be QMA-complete. The Local Consistency problem (deciding 
whether descriptions of small pieces of a quantum system are consistent) is also 
known to be QMA-complete. Here we consider special cases of Local Hamiltonian, 
for "stoquastic" and 1-dimensional systems, that seem to be strictly easier than 
QMA. We show that there exist analogous special cases of Local Consistency, that 
have equivalent complexity (up to poly-time oracle reductions). Our main technical 
tool is a new reduction from Local Consistency to Local Hamiltonian, using SDP 
duality. 

1 Introduction 

Local Hamiltonian is the problem of estimating the ground state energy of a quan- 
tum system with local interactions. This is an important problem in condensed matter 
physics and quantum chemistry, and it is an interesting "matrix-valued" generalization 
of classical problems like Max-k-SAT. Local Hamiltonian also plays a significant role 
in complexity theory: it is QMA-complete, where QMA (Quantum Merlin-Arthur) is a 
generalization of NP where one allows the witness to be a quantum state |17[ El [TBI [24] . 

Two special cases of Local Hamiltonian are of particular interest, because they illus- 
trate the differences in complexity between "quantum" and "classical" problems. First, 
consider a 1-D system: a collection of (i-dimensional particles (qudits) arranged on a 
line, with nearest neighbor interactions. Such systems have been extensively studied in 
condensed matter physics, and in many cases they can be solved efficiently using heuris- 
tic methods (see, e.g., [271 125] ). On the other hand, this problem was recently shown to 
be QMA-hard for d > 12 [21 [J3], and its complexity for smaller d is still open. This com- 
plicated picture is very different from the classical situation: when restricted to a 1-D 
chain, nearly all classical constraint satisfaction problems can be solved exactly in poly- 
nomial time, using dynamic programming. The reason seems to be that 1-D quantum 
systems can have much a richer correlation structure, due to quantum entanglement, 
than 1-D classical systems. 

Secondly, consider a "stoquastic" quantum system [Hi [8]. Here, the Hamiltonian has 
a certain generic form, which ensures that the ground state is a superposition of classical 
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basis states with real non-negative coefficients; thus it more closely resembles a classical 
probability distribution. (This feature is found in many physical systems, as well as 
some quantum algorithms |10|.) Stoquastic Local Hamiltonian is in AM, which suggests 
that it is strictly easier than QMA (unless QMA is in AM, which would be somewhat 
surprising); however, it is also MA-hard, hence it is at least as hard as classical NP- 
complete problems [9]. Also, Stoquastic Local Hamiltonian is complete for a peculiar 
class, StoqMA [8]. Thus stoquastic systems occupy a curious middle ground between 
the quantum and classical regimes. 

In this paper we approach these questions from a different direction. We consider the 
Local Consistency problem, also called consistency of density matrices [T] : given several 
descriptions of small pieces of a quantum system, decide whether they are consistent 
with a single overall state. This problem is QMA-complete under poly-time oracle 
reductions |18j . (A related problem in quantum chemistry, 7V-representability, is also 
QMA-complete [20].) Thus, the Local Consistency and Local Hamiltonian problems 
have equivalent complexity, up to poly-time oracle reductions. In this paper we show 
that this equivalence also holds for special cases involving 1-D and stoquastic systems, 
which are not known to be QMA-hard. 

Our results are as follows. We define special versions of the Local Consistency prob- 
lem for 1-D and stoquastic systems, and, in these special cases, we give poly-time oracle 
reductions from Local Hamiltonian to Local Consistency and vice versa. The reductions 
in the forward direction are similar to [IB] , but the backwards reductions use a new idea 
based on semidefinite programming duality. (Note that, since these problems are not 
known to be QMA-hard, we cannot use the machinery of QMA-completeness; instead 
we have to show reductions in both directions explicitly.) 

This duality idea is interesting in its own right, as it shows a basic connection between 
the Local Hamiltonian and Local Consistency problems. Specifically, we show a poly- 
time oracle reduction from Local Consistency to Local Hamiltonian, that preserves the 
structure of the quantum system (i.e., the number of qubits and the subsets of qubits 
which interact). As mentioned earlier, this reduction also works in cases that are not 
QMA-hard. 

This is quite different from the reduction that is obtained using QMA-completeness. 
There, one gets a poly-time mapping reduction, from Local Consistency to a version of 
Local Hamiltonian that is QMA-hard. That reduction does not preserve the structure of 
the quantum system; the reason is that the Hamiltonian works by "simulating" the QMA 
verifier for Local Consistency, which requires many additional qubits. (In particular, the 
witness for Local Consistency contains multiple copies of the original system [18] , and 
additional qubits are used to simulate the verifier's computation [T71 l3j.) 

Independently, a similar duality technique was used in [12] to study a related problem, 
known as "subsystem compatibility." There, one is given descriptions of large subsets 
of the system, i.e., all but one of the qubits. For a system of n qubits, the input is 
exponentially large in n, and one can solve the problem in time polynomial in the length 
of the input. In our problem, Local Consistency, one is given descriptions of small 
subsets, i.e., k qubits, where k = 0(1). Then the input is only polynomially large in n, 
and in polynomial time, one gets a reduction to Local Hamiltonian instead. 

Outline of the paper: In section 2 we review the definitions of the problems and the 
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basic tools used in the reductions. In section 3 we present our new reduction from Local 
Consistency to Local Hamiltonian via SDP duality. In section 4 we apply this to 1-D 
systems, and in section 5 we apply this to stoquastic systems. 

2 Preliminaries 

2.1 Definitions of Problems 

We briefly review some basic definitions regarding quantum states; see |23j for a thorough 
introduction. Suppose we are interested in a system of n qubits. This is associated with 
a vector space 7i = C 2 ™. 7i has a tensor product structure, 7i = C 2 <g> • • • <g> C 2 (n 
times), and an orthonormal basis consisting of the vectors \z) = \z\) ® ■ ■ ■ <g) \z n ), where 
z\, . . . , z n G {0, 1} (these correspond to the classical states of the system). 

An operator on the system is a matrix acting on Ti. An operator on a subset of qubits 
C C {1, . . . , n} is a matrix of the form M <g> I, where M acts on the tensor factors of TC 
corresponding to C, and / acts on the tensor factors corresponding to {1, . . . , n} — C. 

A quantum state is represented by a density matrix p, which is a positive semidefinite 
matrix with trace 1 acting on TC. A subset of qubits C is described by a reduced density 
matrix tr^ n i_ c (p) (called the partial trace of /?); this is analogous to computing a 
marginal probability distribution by integrating over the unwanted variables. 

The Local Hamiltonian problem is defined as follows \17\ [3] : 

Consider a system of n qubits. We are given a Hamiltonian H = H\ + • • • + 
H m , where each Hi acts on a subset of qubits Ci C {1, . . . ,n} (and so has 
dimension 2^* x 21*^1). The Hi are Hermitian matrices, with norm \\Hi\\ < 1. 
Also, each subset Cj has size |Cj| < k, for some fixed k. 

In addition, we are given a string "I s " (the unary encoding of a natural 
number s), and two real numbers a and 6, such that b — a > 1/s. 

All numbers are specified with poly(m, s) bits of precision. 

The problem is to distinguish between the following two cases: 

• If H has an eigenvalue that is < a, output "YES." 

• If all the eigenvalues of H are > b, output "NO." 

Note that the Hamiltonian H may contain multiple terms that act on the same 
subset, i.e., the subsets Cj might not all be distinct. 

The string "I s " is simply a device to ensure that the gap between the "YES" and 
"NO" cases is not too small, relative to the "size" of the problem. Intuitively, we 
are interested in instances where A: is a constant, m < poly(n) and s < poly(n) (so 
b — a > l/poly(n)). Then we can think of n as the size of the problem, and we say an 
algorithm is efficient if it takes time poly(n). 

Some special cases of this problem are 2-Local Hamiltonian (where k = 2), and 2- 
Local Hamiltonian on a graph G (where k = 2, and the graph G' , consisting of vertices 
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1, . . . , n and edges Ci, . . . , C m , is restricted to be a subgraph of G). In many cases the 
problem has been shown to be QMA-hard [III [111 [2ll El [22l HS1 [26] . 

We remark that classical problems such as Max-k-SAT correspond to the special case 
of Local Hamiltonian where each Hi is a diagonal matrix. 

We define the Local Consistency problem as follows pQ: 

Consider a system of n qubits. We are given a collection of local density 
matrices pi, . . . ,p m , where each pi acts on a subset of qubits Ci C {1, . . . , n} 
(and so has dimension 2^1 x 2'°*'). Each subset Ci has size |Cj| < k, for 
some constant k. 

In addition, we are given a string "I s " (the unary encoding of a natural 
number s), and a real number 0, such that > 

All numbers are specified with poly(s) bits of precision. 

The problem is to distinguish between the following two cases: 

• There exists an n-qubit state a such that, for all i, tr.n ... >n }-Ci (°0 = Pi- 
In this case, output "YES." 

• For all n-qubit states a, there exists some i such that ||tr.r 1; ,.. )n }_c? j (c) — 
Pi 111 ^ in this case, output "NO." 

Without loss of generality, we can assume that the subsets Ci are all distinct; thus 
m < (^) <n k . As in the Local Hamiltonian problem, the string "I s " is simply a device 
to ensure that the gap between the "YES" and "NO" cases is not too small, relative to 
the "size" of the problem. Here, we use the norm \\A\\\ = tr \A\ to measure the distance 
between pi and the corresponding reduced density matrix of a. When multipled by 1/2, 
this is the trace distance. 

An important special case is where k = 2, and we can visualize the system as a graph 
with nodes 1, . . . , n and edges given by the subsets C±, . . . , C m . Local Consistency was 
shown to be QMA-hard in [181 I19j . via an oracle reduction from Local Hamiltonian. 
This reduction preserves the subsets C±, . . . , C m ; thus, for many of the special cases in 
which Local Hamiltonian is QMA-hard, Local Consistency is also QMA-hard. 

We remark that when the matrices pi are all diagonal, this reduces to a classical 
problem of deciding the consistency of marginal probability distributions. 

In sections 4 and 5 we will define special cases of Local Hamiltonian and Local 
Consistency, for 1-D and stoquastic systems. 

2.2 Convex Optimization using a Membership Oracle 

We review some algorithms for convex optimization using a membership oracle |28} ITT1 
114] . which were the main tool in reducing Local Hamiltonian to Local Consistency 
[18[ [19] . This is a summary of the results in [19] , which were based on |llj . 

First, some notation: let S(p, r) denote the closed ball of radius r around the point 
p, S(p, r) = {x G M. n | ||x — p\\ < r}. Also, for any set K, we define the ball of radius e 
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around K, 

S{K,e) = {x G R n | there exists y G K s.t. ||x - y\\ < e}, 
and we define the interior of K with radius e, 

-e) = {x G M n | S(x, e) C K}. 

Let X be a closed convex set in M n , and suppose we are given a point p G M n , and 
inner and outer radii r,iJ€l, such that S(p,r) C if C S(0,R). (This implies that 
is bounded and full-dimensional.) We want to show a reduction from the problem of 
optimizing a linear function over K, to the problem of deciding membership in K. 

We define the weak optimization problem WOPT £ as follows: (The adjective "weak" 
refers to the fact that we allow additive errors of size e.) 

Given c G M n , ||c|| = 1, 7 G R, and e G R, e > 0, all specified with poly(n) 
bits of precision. 

If there exists a vector y G S(K, —e) with c • y > 7 + e, then answer "YES." 
If for all x G e), c • a? < 7 — e, then answer "NO." 

We define the weak membership problem WMEMg as follows: 

Given y G W 1 , and 5 G R, <5 > 0, all specified with poly(ra) bits of precision. 
If y G S(K, -5), then answer "YES." 
If y <£ S(K,6), then answer "NO." 

A reduction from WOPT £ to WMEM$ is given in [11], using the shallow-cut ellipsoid 
method. However, we need to modify their result slightly. The reduction in [11] uses 
a model for "exact" convex optimization, where e is exponentially small, and 5 may be 
exponentially smaller than e. But for our application, e is inverse-polynomial in size, 
and we want 5 to be at most polynomially smaller than e; this can be described as 
"approximate" convex optimization. 

It turns out that the reduction of [11] will suffice for our purposes, as long as R/r is 
at most polynomially large. The following statement is proved in |19j . 

Theorem 1 Let K be any closed convex set in R n , such that S(p,r) C K C S(0,R), as 
defined above. Then there is an oracle reduction from WOPT £ to WMEM$, for some 
5 > poly(e, (r/R), (1/n)), which runs in time poly(n, (R/r), (1/e)). 

We state a straightforward corollary of this result, that will be more convenient for 
our purposes. First, we define a slightly modified problem, WOPT*, as follows: 

Given c G M n , ||c|| = 1, 7 G R, and e G R, e > 0, all specified with poly(n) 
bits of precision. 

If there exists a vector y G K with c • y > 7 + e, then answer "YES." 
If for all x G K, c ■ x < 7 — e, then answer "NO." 

We also define a new problem, WMEM$ , as follows: 
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Given y £ W 1 , and 5 € M, 5 > 0, all specified with poly(n) bits of precision. 
UyeK, then answer "YES." 
If y <£ S(K,5), then answer "NO." 

Corollary 2 Let K be any closed convex set in W n , such that S(p, r) C K C 5(0, R), as 
defined above. Then there is an oracle reduction from WOPT* to WMEM^, for some 
5 > poly(e, (r/R), (1/n)), which runs in time poly(ra, (R/r), (1/e)). 

2.3 Reduction from Local Hamiltonian to Local Consistency 

We briefly review the main result of |18| . 

Theorem 3 There is a poly-time oracle reduction from Local Hamiltonian to Local Con- 
sistency. 

The main idea is that finding the lowest eigenvalue of a local Hamiltonian H = 
H\ + • • • + H m is equivalent to finding local density matrices p%, . . . , p m that minimize 
the linear function /(pi,...,p m ) = tv(Hipi) + ••• + tv(H m p m ) over the convex set 
K = {(pi, . . . , p m ) that are consistent}. (Note that, if pi, . . . ,p m are consistent with a 
global state a, then /(pi, . . . , p m ) = tv(Ha).) Given an oracle for the Local Consistency 
problem, we can construct a membership oracle for the set K, then apply Corollary [2] to 
solve the optimization problem over K, and thus solve the Local Hamiltonian problem. 

The main technical detail is to formulate the problem in such a way that K has 
the necessary geometric properties, i.e., K has outer radius R and inner radius r such 
that R/r is at most polynomially large. To this end, we will construct a set 5 of local 
observables; their expectation values encode the information contained in pi,...,p m 
without any redundancy. 

We introduce the n-qubit Pauli matrices P = ®" =1 Pi, where Pi £ {I, X, Y, Z}. Note 
that for any n-qubit Pauli matrices P and Q, tr(PQ) = 2 n if P = Q and otherwise. 
Any n-qubit density matrix a can be written in the form 

° = ^Z^2 apP > a P = tr(Pa). 
p 

Let C be a subset of qubits. We say that P is supported inside C if, for all i ^ C, 
Pi = I. We claim that the expectation values ap, for those P supported in C, contain 
precisely the same information as the reduced density matrix on C. To see this, define 
P\c = t^igc^' which we call the "restriction" of P to C. Then we can write the 
reduced density matrix on C in the form 

tr{i,...,»}-c(<r) = ^ Yj a P tT {i,...,n}-c( p ) = 7^c\ Yl a P p \c- 

P supp. in C P supp. in C 

For a collection of subsets C\,... ,C m , we define 5 to be the set of "local" Pauli 
matrices, excluding the identity matrix, 

m 

S = \^J{P | P is supported inside Ci} — {I}. 

i=l 
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We also let D = \S\, and note that D < A k m — 1. We can now replace the local density 
matrices pi, ■ ■ ■ ,p m with the set of expectation values otp for the observables P G S. 
Using this formulation, the convex set K has the required geometric properties. (This 
uses the orthogonality properties of P; see [18] for details.) 

3 Reduction from Local Consistency to Local Hamiltonian 

The idea comes from a theorem of "strong alternatives" in semidefinite programming [5] . 
Let F\, . . . ,Fd be complex Hermitian matrices of dimension N. Consider the following 
matrix inequality: 

D 

Y^XiFi + I < 0, (1) 

i=l 

where x G ~R. D is a variable. (Notation: M -< means M is strictly negative definite, 
M >z means M is positive semidefinite, etc.) Also consider the following system of 
inequalities: 

Z±0, Z^O, tv(FiZ) = 0(Vi = l,...,D), (2) 

where Z, a complex Hermitian matrix of dimension N, is a variable. The theorem states 
that exactly one of the two inequalities (pQ) and ([2]) is feasible. In other words, if ([2]) 
is feasible, then (pQ) is not; and if ([2]) is not feasible, then ([I]) is. (When this property 
holds, we say that (H|) and ([2]) are strong alternatives.) 

Observe that inequality ([2]) can be used to express the Local Consistency problem: 
Z is a global density matrix (unnormalized, but note that all the constraints remain the 
same if we divide across by tr(Z)), and we can choose the constraints ti(FiZ) = to 
ensure that Z agrees with the desired local density matrices (note that the matrices F{ 
will then be local observables). But now the expression ^2d = i x i^i + i n inequality ([T]) is 
simply a local Hamiltonian, and estimating its largest eigenvalue is precisely the Local 
Hamiltonian problem (modulo a sign flip). So a Local Hamiltonian oracle allows us to 
test membership in the convex set defined by inequality ([T]); and, using the methods 
of convex optimization described in Chapter 2, we can then decide the feasibility of 
([T]). Since ([T]) and ([2]) are strong alternatives, this lets us solve the Local Consistency 
problem. 

This is the intuition, but some further work is needed to make it rigorous. We 
have to allow for the inverse-polynomial precision in the Local Consistency and Local 
Hamiltonian problems. Also, in order to do convex optimization with a membership 
oracle, the set of feasible solutions K must satisfy certain geometric properties. (In 
particular, K must be finite!) So we have to formulate inequality JI| in a different 
way. We will show a finite-precision, "algorithmic" version of the theorem of strong 
alternatives. 

Theorem 4 There is a poly-time oracle reduction from Local Consistency to Local 
Hamiltonian. 

Proof: Suppose we have an instance of the Local Consistency problem, i.e., a collection 
of local density matrices p%,...,p m , where pi describes a subset of qubits Cj, and an 
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error parameter (5. Let S be the set of local Pauli observables, and let D = \S\. For each 
observable P € S, we define ap to be the desired expectation value, which we compute 
as follows: pick some subset Cf such that P is supported in Cj, then set ap = tr(Ppi). 

We write down a convex program and its dual. For each P £ S, we define a new 
observable 

F P = P-a P I, 

which is shifted so that the desired expectation value now equals 0. We also define F{x) 
to be a linear combination of these observables, 



F(x) = ^ XpFp + I, for x £ 



!>D 

h pi: p t i , iui x c 14 . 

Pes 

Now consider the following convex program: 



Find some x £ [-1, 1] D and s € [1 - 2D, 1 + 2D] that 
minimize s such that F{x) < si. 

To see that this is a convex program, recall that the largest eigenvalue of F(x) is a 
convex function of x, since it can be written as the pointwise minimum over a family of 
affine functions of x. The variable s is redundant here, but it will play a role later when 
we apply algorithms to solve this program. We will refer to this as the primal program; 
let p* denote the optimal value of the objective function s. 
The dual program is as follows: 

Find some 2 n x 2 n complex matrix Z that 
maximizes g(Z) such that Z >z and tr(Z) = 1, 

where the dual function g(Z) is given by 

g(Z) = inf 8 + tv(Z(F(x) - si)) 

xe[-i,i] D 

se[l-2D,l+2D] 

= inf tr(ZF(x)) 
xe[-i,i] D 

= inf V xptr(ZFp) + 1. 

Let al* denote the optimal value of the objective function g{Z). Strong duality holds 
because the primal problem is convex and satisfies a generalized Slater condition [5] (to 
see this, note that the point (x,s) = (0,2) is strictly feasible, i.e., it lies in the relative 
interior of the domain, and it satisfies F(x) -< si). Strong duality implies that p* = d* , 
i.e., the optimal values of the primal and dual programs are equal. 

We now give a poly-time oracle reduction from Local Consistency to Local Hamil- 
tonian. We show that Local Consistency reduces to the weak optimization problem 
WOPT*, which reduces to the weak membership problem W MEM* , which reduces to 
Local Hamiltonian. 

First, suppose we have a "YES" instance of Local Consistency. Then there exists an 
n-qubit state a such that, for all P £ S, tr(Pcr) = ap. So in the dual program there 
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exists some Z y 0, tr(Z) = 1, such that for all P £ S, tr(ZFp) = 0. This implies 
g(Z) = 1, hence the dual program has optimal value d* > 1. By strong duality, the 
primal program has optimal value p* > 1. 

On the other hand, suppose we have a "NO" instance of Local Consistency. We 
claim that, for all n-qubit states a, J^PeS I tr(Pa) — ap| > /?. This can be seen as 
follows. Note that, for any a, there is some subset C{ such that \\a — pi\\\ > (3, where 
a = tip Using the matrix Cauchy-Schwarz inequality [BJ, ||cr — Pi||i < 

|| cr — pj[|2"v2*. By Fourier analysis, 

P supp. on Ci 

s4[EifW*-«»i 

vz PeS 

The claim follows by combining these inequalities. 

Therefore, in the dual program, for all Z such that Z >z and tr(Z) = 1, we have 
that Xlpgs I tr(-^-Pp)| > A which implies g(Z) < 1 — /?. Thus the dual program has 
optimal value d* < 1 — /?. By strong duality, the primal program has optimal value 
p* < 1 - p. 

So we have reduced Local Consistency to the problem of distinguishing between the 
two cases p* > 1 and p* < 1 — (5 for the primal program. This is an instance of the weak 
optimization problem WOPT^, 2 over the convex set 

K = {(a?, s) £ [-1, 1] D X [1 - 2D, 1 + 2L>] | F(sc) - si ^ 0}. 



1/2 



A)) 2 

-=^|tr(Pa)-a P |. 
vz Pe<s 



Now we will reduce WOPT* to W MEM* . First we need some bounds on the 
geometry of K. It is easy to see that K is contained within a ball of radius R = 
wD + (1 + 2D) 2 < 0(D). In addition, we claim that K contains a ball around the 
point (0, . . . , 0, 2) of radius r = 4 ^ 1 +1 ^ . To see this, consider an arbitrary point (y, t + 2) 

where y £ R D , t £ R and y^yll 2 + t 2 < 4(D 1 +1) . The operator 

F{y) - (t + 2)1 = £ ypFp — tl — I 
Pes 

has all of its eigenvalues bounded above by X^pgs 4(j-fi) II-^-pII + A(D+i) ~ ^ — ~ 2 ( usm g 
the fact that ||-Fp|| < 2). Thus (y,t + 2) is in K. This proves the claim. 

So we have 22/r < 0(D 2 ). By Corollary El WOPT* /2 reduces to WMEM* S where 
5 > poly(/3, 1/D), with running time poly(D, 

Finally, we reduce W MEM* to the Local Hamiltonian problem. Observe that, since 
the Fp are local operators, F(x) is a local Hamiltonian. Given an oracle that solves 
the Local Hamiltonian problem, we can estimate the largest eigenvalue of F{x) (i.e., the 
smallest eigenvalue of —F(x)), and thus decide whether (x,s) is in the set K. 

Suppose we have a "YES" instance of WMEM* & . Then (x,s) £ K, so F(x) < si, 
i.e., all eigenvalues of —F(x) are > —s. So this is a "NO" instance of Local Hamiltonian. 
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Now suppose we have a "NO" instance of W MEM$ . Then (x,s) ^ S(K,8), and in 
particular, (x, s + 5) ^ K. So F(x) ^ (s + 8)1, i.e., —F(x) has an eigenvalue that is 
< —s — 5. So this is a "YES" instance of Local Hamiltonian. 

Note that \\F(x)\\ < Y.p e sW F P\\ + 1 < 2D + 1. Thus, WMEM* S reduces to Local 
Hamiltonian with precision 5/ {2D + 1). 

Thus we conclude that Local Consistency (with precision 0) reduces to Local Hamil- 
tonian (with precision poly(/3, 1/D)), and the running time is poly(-D, 1//?). Note that 
D < A k m is polynomial in the size of the input. □ 

4 Local Consistency for 1-D Systems 

Let us consider a 1-dimensional chain of n qudits (a qudit is a d-dimensional particle), 
with nearest-neighbor interactions (i.e., interactions between particle i and particle i + 1, 
for i = 1, . . . , n — 1). 

First consider the case of qubits {d = 2). The reduction from Local Hamiltonian to 
Local Consistency shown in Theorem El and the reverse reduction shown in Theorem HJ 
both preserve the neighborhood structure of the problems — that is, each local term in 
the Hamiltonian corresponds to a local density matrix, and vice versa. Thus we have: 

Corollary 5 On a 1-D chain of qubits (d = 2), Local Hamiltonian and Local Consis- 
tency have equivalent complexity (up to poly-time oracle reductions). 

We will now sketch one way of extending these results to the case of qudits (d > 2). 
The first step is to define a set of observables for a single qudit, with nice properties 
similar to the Pauli matrices. Let \i), i = 0, . . . , d — 1 denote the standard basis states 
for a single qudit. Also, let i (in plain, not italic type) denote the square root of —1. 



Note that Z{ is the diagonal matrix whose diagonal consists of in the first i + 1 
positions, followed by —1, followed by in all the remaining positions. We have a total 
of 2g) + (d - 1) = d(d - 1) + (d - 1) = d 2 - 1 observables. 

These observables satisfy the following orthogonality relations: 




0<i<j <d-l 
<i < j < d-1 




0<i<d-2 
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tr(AB) 
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I 


Xki 





I 


Y kl 





I 


Zk 







Xki 


2 if = (k,l); otherwise 


x^ 


Yki 





x^ 


Zk 





Y 


Yki 


2 if = (k,l); otherwise 


Y 


Zk 





Zi 


Zk 


1 + t-^-t if i = k; otherwise 



In addition, note that = = \\Zi\\ = 1. 

We can now use these qudit observables in the same way that we used the Pauli 
matrices for qubits. We construct n-qudit observables by taking tensor products of 
single-qudit observables: P = (^)^ =1 P a , where P a G {I, Xij,Yij, Zi}. Note that for any 
n-qudit observables P and Q, tv(PQ) = tr(P 2 ) if P = Q and otherwise. Any n-qudit 
density matrix a can be written in the form 



x -« ap 



^tr(J«) 



P, a P = tr(Pcr). 



As before, we say that P is supported inside a subset C C {1, . . . , n} if for all i ^ C, 



Pi = I. We define P\ c 



Pi, the "restriction" of P to the subset C . We can write 



the reduced density matrix for the subset C in the form 

ap 



tr 



{l,...,n}-c(< 7 ) 



E 



P supp. in C 



tr(P 2 ) 



tr 



{l,...,n}-C 



(P) 



E 



ap 



P supp. in C 



tr((P|c) 2 )' 



\C- 



Now we can use essentially the same reductions as before, from Local Hamiltonian to 
Local Consistency and vice versa, for systems of qudits. (Details omitted.) In particular, 
this implies: 

Corollary 6 On a 1-D chain of qudits (for any fixed d > 2), Local Hamiltonian and 
Local Consistency have equivalent complexity (up to poly-time oracle reductions). 



5 Stoquastic Local Consistency 

A Hamiltonian H is called "stoquastic" if all of its off-diagonal matrix elements, relative 
to the standard basis, are less than or equal to 0. (Note that the diagonal elements can 
be made < by subtracting a multiple of the identity from H; this shifts the eigenvalues 
but does not change the eigenvectors.) By the Perron- Frobenius theorem [4], the ground 
state of H has the form = ^2 Z c z \z) (up to an overall phase factor), where \z) are the 
standard basis vectors and c z > 0. (In other words, \ip) is a superposition without any 
negative or complex coefficients; thus it resembles a classical probability distribution.) 

In the Stoquastic Local Hamiltonian problem, we are given a local Hamiltonian 
H = H\ + • • • + H m where each of the local terms Hi is stoquastic. As discussed 



11 



previously, this makes the problem potentially easier [9j [8] . In this section we propose 
a Stoquastic Local Consistency problem, and show that it has the same complexity as 
Stoquastic Local Hamiltonian (up to poly-time reductions). 

First, let us say that a density matrix p is "stoquastic" if all of its off-diagonal matrix 
elements, relative to the standard basis, are greater than or equal to 0. (Its diagonal 
elements must be > since p is positive semidefinite.) Note that the set of stoquastic 
density matrices is convex. 

Now consider an obvious way of defining the Stoquastic Local Consistency problem: 

Given local density matrices pi, . . . , p m which are stoquastic, does there exist 
a global density matrix a that agrees with pi, . . . , p m ? 

Stoquastic Local Hamiltonian easily reduces to this problem. However, it is not clear 
whether this problem in turn reduces to Stoquastic Local Hamiltonian; the technique 
from section [3] instead produces a reduction from this problem to standard Local Hamil- 
tonian. 

Instead we will use a more subtle definition of Stoquastic Local Consistency: 

Given local density matrices pi,---,p m , does there exist a global density 
matrix a such that, for all i = 1, . . . , m, tr.n ... %n }-Ci ( CT ) — e Pi? 

(Here, Ci is the subset of qubits described by pi, and > e denotes element-wise 
inequality between two matrices written in the standard basis; we assume all 
matrices are real.) 

This definition is a little unusual, but it turns out to have the desired property: we can 
give a reduction from Stoquastic Local Hamiltonian to Stoquastic Consistency, and vice 
versa. The element-wise inequality > e comes about naturally from the duality technique 
used in the second reduction. 

We now state the full definition of the Stoquastic Local Consistency problem: 

Consider a system of n qubits. We are given a collection of local density 
matrices pi, . . . , p m , where each pi acts on a subset of qubits Ci C {1, . . . , n}. 
The matrices are assumed to be real (not complex). Also, each subset Ci has 
size \Ci\ < k, for some constant k. 

In addition, we are given a string "I s " (the unary encoding of a natural 
number s), and a real number f3, such that (3 > 1/s. 

The problem is to distinguish between the following two cases: 

• There exists a real n-qubit state a such that, for all i, and for all s, t G 

{0, l}' Cl ', we have 

(s|tr { i,...,n}-C l O)|£} > (s\pi\t)- 
In this case, answer "YES." 

• For all real n-qubit states a, there exists some i, and there exist some 
s,t G {0, ljl^l, such that 

In this case, answer "NO." 
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We will show the following result: 

Theorem 7 Stoquastic Local Hamiltonian and Stoquastic Local Consistency have equiv- 
alent complexity (up to poly-time oracle reductions). 

5.1 Prom Stoquastic Local Hamiltonian to Stoquastic Local Consis- 
tency 

First we show a reduction from Stoquastic Local Hamiltonian to Stoquastic Local Con- 
sistency. The basic idea is as follows. We are given a local Hamiltonian H = 
where the Hi are real and stoquastic. Without loss of generality, we can assume Hi < e 
(we simply add a multiple of the identity to Hi). Also, let C L be the subset of qubits on 
which Hi acts. 

Let px, . . . ,p m be local density matrices, where pi describes the subset of qubits Cj. 
We want to find p%, . . . , p m that correspond to the ground state of H. We consider the 
following convex program: 

Find pi, . . . , p m that minimize YliLi tr(-HiPi), subject to the constraints: 

1. For all i, pi >z and tr(/>j) = 1. 

2. There exists an n-qubit state a s.t. a >z 0, tr(cr) = 1, 
and for all i, tr {1 . ^.^(ct) > e p». 

Here, all matrices are restricted to be real. 

We claim that this convex program is equivalent to the Stoquastic Local Hamiltonian 
problem. If H has an eigenstate \ip) with eigenvalue < A, then the convex program has 
optimal value < A; to see this, set p. L = n }-Ci IvK^I- O n the other hand, if all the 

eigenvalues of H are > A + 5, then the convex program has optimal value > A + 5; this 
follows because, for any feasible pi, . . . , p m , we have that 

m m 

J2^(H iPi ) > J^tr(fli<r) = tr(Ha), 

i=l i=i 

using constraint (2) and the fact that Hi < e 0. 

Solving this convex program is an instance of the optimization problem WOPT* 
over the convex set K defined by constraints (1) and (2). By Corollary [21 this reduces 
to the membership problem W MEM* over K. This immediately reduces to Stoquastic 
Local Consistency. 

The main technical detail is to formulate the problem so that the set K is full- 
dimensional, with inner and outer radii that satisfy R/r < poly(n). This can be done 
using a subset of the local Pauli observables, as in Theorem [31 with two minor modifi- 
cations. 

First, in this problem the matrices p±, . . . , p m are restricted to be real. To represent 
these, we use the subset of local Pauli observables P = ®ILi ^ that contain an even 
number of Y factors. This follows from the equation 

a= 2n^2 apP ' a P = tv(Pa); 
p 
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note that ap is always real (since a is Hermitian), while P is real (imaginary) when it 
contains an even (odd) number of Y factors. One can check that there are (4 n + 2 n )/2 
real observables P, which correspond to the degrees of freedom for a 2" x 2™ real matrix. 

Secondly, in this problem the set K is larger than in Theorem O since the constraint 
(2) is weaker. Nonetheless, constraint (1) ensures that K is still contained within a ball 
of radius R. K also contains a ball of radius r, by the same argument as before. 

5.2 Prom Stoquastic Local Consistency to Stoquastic Local Hamilto- 
nian 

Next we show a reduction from Stoquastic Local Consistency to Stoquastic Local Hamil- 
tonian. The reduction uses strong duality, as in Section [3j 

Suppose we have an instance of Stoquastic Local Consistency. The first step is to 
represent pi, ... ,p m as the expectation values of certain observables. However, we use 
a different set of observables, instead of the Pauli matrices, so that we can deal with 
inequalities involving the matrix elements of pi. These observables do not have any nice 
orthogonality properties, but the reduction technique from Section [3] does not require 
this. 

For each subset of qubits Cj, define the following observables acting on Cf 
X$ = ±(\s)(t\ + \t)(s\), fl ,t€{0,l}l<*l, s±t, 

where s ^ t denotes lexicographic order. We can think of these observables as acting on 
the full n-qubit system (we tensor them with the identity matrix). For any real n-qubit 
state a, the matrix elements of ton n \_(j.(cr) are given by the expectation values of 
these observables: 

trp^jV) = (s\ti {1 _ n} _ Ci {a)\t). 

Then the conditions for a "YES" instance of Stoquastic Local Consistency can be written 
as: 

tr(x2V) > (s\pi\t). 

(i) 

We let S be the set of all these observables X st , for all of the subsets Cj, i = 1, . . . ,m. 
We also let D = \S\. Note that < 1. 

Next, we formulate a convex program, together with its dual. Define new observables 

F®=X®-(8\pi\t)I, 

which are shifted so that our goal is to satisfy the inequalities 

tr(F«a) > 0. 

For notational convenience, let us refer to these observables as F p , for p = 1, ...,£?. 
Define F(x) to be a linear combination of these observables, 

D 

F(x) = x p F p + I, for x G R D . 
p=i 
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We construct a convex program which is similar to the one in Section [3j except that 
we restrict x to lie in the domain [0,1] D instead of [— 1,1] D . (This restriction of the 
domain is the key feature that will eventually connect the dual program to Stoquastic 
Local Consistency.) 

Find some x G [0, 1] D and s G [1 - 2D, 1 + 2D] that 
minimize s such that F(x) ^ si. 

This is the primal program; let p* denote the optimal value of the objective function s. 
The dual program is as follows: 

Find some 2 n x 2 n real matrix Z that 
maximizes g(Z) such that Z >z and tr(Z) = 1, 

where the dual function g(Z) is given by 

D 

g(Z)= inf tv(ZF(x))= inf V x p tr(ZF p ) + 1. 
*e[o,ip ^[o,i] D ^ 

Let d* denote the optimal value of the objective function g(Z). Strong duality holds 
because the primal problem is convex and satisfies a generalized Slater condition [5] (to 
see this, note that the point (x,s) = ((1/3D)1,2) is strictly feasible). Strong duality 
implies that p* = d*. 

Now, suppose we have a "YES" instance of Stoquastic Local Consistency. Then in 
the dual program there exists some Z >z 0, tr(Z) = 1, such that for all p, ti(ZF p ) > 0. 
This implies g(Z) = 1, hence the dual program has optimal value d* > 1. By strong 
duality, the primal program has optimal value p* > 1. 

On the other hand, suppose we have a "NO" instance of Stoquastic Local Con- 
sistency. Then for all Z such that Z y and tr(Z) = 1, there is some p such that 
tr(ZFp) < —f3, which implies g(Z) < 1 — /3. Thus the dual program has optimal value 
d* < 1 — p. By strong duality, the primal program has optimal value p* < 1 — (3. 

Thus we have reduced Stoquastic Local Consistency to the primal problem. This is 
an instance of the optimization problem WOPT* over the set 

K = {(x, s) G [0, 1] D x [1 - 2D, 1 + 2D] such that F(x) * si}. 

Using a similar analysis as in section [31 we can show that K has inner and outer radii 
that satisfy R/r < poly(D). (For instance, one can center the inner ball around the point 
((1/3.0)1,2).) By Corollary [21 this problem then reduces to the membership problem 
WMEM* for K. 

We claim that we can solve this problem, given an oracle for Stoquastic Local Hamil- 
tonian. Observe that the F p are local operators, whose off-diagonal elements are all > 0. 
Thus —F(x) is a stoquastic local Hamiltonian, and we can use the oracle to estimate 
its ground state energy. This is equivalent to estimating the largest eigenvalue of F(x), 
which allows us to test whether the constraint F(x) ^ si is satisfied, and thus test 
membership in K. 
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6 Discussion 



We have shown that Local Hamiltonian and Local Consistency have equivalent com- 
plexity in two special cases, 1-D chains of qudits with small d and stoquastic quantum 
systems, where neither problem is known to be QMA-hard. To do this, we used ideas 
from |18j . together with a new reduction from Local Consistency to Local Hamiltonian 
using semidefinite programming duality. 

In practice, one would like to solve particular instances of Local Hamiltonian that 
arise in condensed matter physics or quantum chemistry. Reducing the problem to Local 
Consistency is one possible approach. (In fact, this is the underlying idea in 2-RDM 
methods in quantum chemistry [21] .) However, since Local Consistency is QMA-hard, 
it might seem that this only leads to a more general (hence potentially harder) problem. 
Our results show that this is not always the case: for 1-D and stoquastic systems (without 
any other special features), an efficient solution to Local Hamiltonian implies an efficient 
solution to Local Consistency. It would be interesting to find more instances where Local 
Consistency is a potentially tractable problem. 

Finally, one might ask whether Stoquastic Local Consistency sheds any light on the 
class StoqMA introduced in [Sj. (Recall that Stoquastic Local Hamiltonian is StoqMA- 
complete [8].) Clearly, Stoquastic Local Consistency is StoqMA-hard, but we were not 
able to show that it is in StoqMA. Could there be another complexity class that better 
describes these problems? 

Acknowledgements: Thanks to Frank Verstraete and Daniel Nagaj for useful discussions. 
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